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j_| ■ Magnetars^'^ are a special class of slowly rotating (period ~5-12 s) neutron stars with ex- 
tremely strong magnetic fields (>10^^ G) - at least an order of magnitude larger than those 
of the "normal" radio pulsars. The potential evolutionary links and differences between these 
two types of objects are still unknown; recent studies, however, have provided circumstan- 
tial evidence connecting magnetars with very massive progenitor stars Here we report 
the discovery of an infrared elliptical ring or shell surrounding the magnetar SGR 1900+14. 



The appearance and energetics of the ring are difficult to interpret within the framework 
of the progenitor's stellar mass loss or the subsequent evolution of the supernova remnant. 
We suggest instead that a dust-free cavity was produced in the magnetar environment by the 
giant flare emitted by the source in August 1998. Considering the total energy released in 
the flare, the theoretical dust-destruction radius matches well with the observed dimensions 
of the ring. We conclude that SGR 1900+14 is unambiguously associated with a cluster of 
massive stars, thereby solidifying the link between magnetars and massive stars. 

Soft gamma repeaters (SGRs) and anomalous X-ray pulsars (AXPs) are the two main classes 
of objects currently believed to be magnetars. Their characteristic ages, derived from their rota- 
tional properties^, indicate a young population - typically a few thousand years old. Although these 
spin-down ages must be treated with caution, the relative youth of SGRs and AXPs is supported 
by the fact that some AXPs are located at the centres of supernova remnants (SNRs). Associations 
between SNRs and SGRs have also been claimed in the literature, but, unlike the AXPs, the SGRs 
are offset from the centres of their proposed SNRs, increasing the likelihood of spurious align- 
ment. As a result, the validity of all SGR-SNR associations has been questioned^. Both AXPs and 
SGRs, for example the AXP CXOU J164710.2-455216 (ref.5) and SGR 1806-20 (ref.8), have 
been linked to clusters of massive stars. The implied progenitor mass of M > 40 — SOM© supports 
theoretical predictions that very massive stars with sufficient metallicity and the corresponding 
high mass loss can still form neutron stars^, contrary to the standard evolutionary picture that such 
massive stars end their lives as black holes. 
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As part of our systematic study of the circumstellar environments of SGRs and AXPs, we 
observed the position of SGR 1900+14 using all three instruments onboard the NASA Spitzer 
Space Telescope in 2005 and 2007. Surprisingly, our 16 /xm and 24 iim images (Fig. lb, c) reveal 
a prominent ring-like structure that is not detected in our 3.6-8.0 /im observations. A formal ellip- 
tical fit to the ring indicates semi-major and semi-minor axes of angular lenghts ~ 36" and ~ 19", 
respectively, centred at right ascension 19:07:14.32 and declination -1-09:19:20.0 ±1.0". This co- 
incides with the radio position of SGR 1900+14 (19:07:14.33 -1-09:19:20.1 ±0.15"), determined 
from observations of a transient synchrotron source associated with the 1998 giant flare^°. We re- 
examined published and archival data of the field around SGR 1900+14 but found no equivalent 
structure at optical, near-lR, radio or X-ray wavelengths. In particular, the lack of detection in the 
radio spectrum, to a limit of -L332MHZ < 4 x 10^^ df^ erg s~^, and the X-ray spectrum, to a limit 
of L2-iokeV < 2.7 X 10^^ dl^ erg (Chandra data archive, http:cxc.harvard.educda), is important 
in determining the nature of the ring. Here di^ parameterizes the distance to the source, according 
to d = 15 di5 kpc. 

The Spitzer images are dominated by the bright emission from two nearby M5 supergiants 
that mark the centre of a compact cluster of massive stars at a distance of 12-15 kpc^^' Although 
it has been suggested that SGR 1900+14 is associated with this cluster^^, an alternative distance of 
5 kpc has also been suggested, on the basis of the hydrogen column density of its X-ray spectrum^"^. 
In addition, the SGR is offset by 12" from the cluster centre, and the visual extinction correspond- 
ing to the X-ray derived hydrogen column density^^ (Ay = 12.8 ± 0.8) is significantly different 
from that deduced optically for the cluster stars^^ (Ay — 19.2 ± 1). Considering both possibilities, 
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we calculated physical sizes for the semi-major and semi-minor axes of the ring of 0.9 x 0.5 pc, 
if SGR 1900+14 lies at a distance of ~ 5 kpc, or 2.6 x 1.4 pc, if it resides at a cluster distance of 
15 kpc. 

We have constrained the temperature of the material in the ring by first cross-convolving the 
images with the point spread function of the different arrays. We then measured the flux through 
fixed apertures at several positions along the ring. The largest source of uncertainty lies in the 
subtraction of the local background. We derived a temperature of 130 — 150 K for the material 
in the ring using a simple black-body fit, and one of ~ 80 — 120 K using a more realistic dust 
model. We caution that the presence of spectral lines or peculiarities in the underlying mid- 
infrared spectral shape might significantly alter these derived temperatures. For example, Spitzer 
observations have revealed a shell only visible at 24 iim, where the broadband flux can be entirely 
attributed to [OIV] line emission^^. A broad 22 //m continuum feature has also been identified in 
environments associated with supemovae and massive star formation^^. 

The 70 fjLm image shows a bipolar flux distribution along the minor axis of the ellipse (Fig. 
Id). Extended enhancements in the ring are also seen at these positions in the 16 and 24 /xm images. 
Our data do not allow us to determine the underlying three-dimensional geometry of the structure, 
that is, to distinguish between a limb-brightened shell, a true ring morphology, or a bipolar cavity 
with equatorial torus. Similarly, an accurate measurement of the total flux in the ring is very 
difficult because the enveloping diffuse emission (Fig. la) coupled with the contribution from the 
M5 supergiants prohibits a clean separation of the ring emission. Using a narrow elliptical aperture. 
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we measured fluxes of 1.2±0.2 Jy and 0.4±0.1 Jy at 24 //m and 16 //m, respectively, implying ring 
luminosities of v L^, {16 jim) ~ 2 x 10^^ dj^ erg and uLy{2AiJim) ~ 4 x 10^^ erg (where 
V denotes the frequency at which the luminosity is determined). 

To test for ring evolution, we used the publicly available MIPSGAL Spitzer Legacy pro- 
gramme observations of the field that were obtained at an epoch (2005 Oct 03) earlier than was our 
own data set. Creating a difference image from the 24 jira data of the two different epochs, we find 
no discernable change in the size (positional shifts < 0.5") or flux (root mean square < 1 — 2%) 
of the ring emission in the 1.6 yr between the two observations. The stationary nature of the ring 
strongly disfavours one obvious explanation for the creation of the ring - a light echo, due to re- 
processing by the surrounding dust, from the 1998 SGR 1900+14 giant flare. This model predicts 
a typical fractional change in the ring size and brightness between the two epoch observations of 
~ 23%, which is not seen in the data (see Supplementary Information for more details). 

Morphologically, the ring resembles the wind-blown bubbles and shells observed around 
evolved massive stars^^"^*^, such as B supergiants, luminous blue variables and Wolf-Rayet stars. 
The supersonic wind from the star drives a shock into the surrounding medium, sweeping up 
the ambient material into a thin, dense, cool shell. The emission from the swept-up material is 
powered by the luminous central star. Such shells have typical radii of ~ 2 — 10 pc for Wolf-Rayet 
stars^^ and 0.1 — 2.3 pc for luminous blue variables/supergiants^^'^^ , overlapping the size of our 
ring. To guard against a simple chance superposition with an unrelated source, we investigated the 
photometric properties of stars near the centre of the ring. The sources in the immediate vicinity 
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of the SGR position^ ^ are too faint in comparison with the near-infrared luminosities expected for 
the massive stars discussed above^'-^^, and thus cannot be physically associated with the ring. For 
stars with larger offsets from the centre (> 7") we cannot rule out a massive star classification 
based solely on available archival data. However, a situation in which an off-centre star creates an 
elliptical ring with an unrelated magnetar at its exact centre appears rather contrived. 

In the absence of an obvious chance superposition, the close positional coincidence between 
the SGR and the centre of the ring then strongly indicates that the SGR and the ring are phys- 
ically connected. One interpretation is that the ring represents material ejected during the late 
stages of the SGR progenitor's evolution. However, we know that a supernova explosion occurred, 
producing the magnetar, and the resulting shock wave is expected to disrupt such a close shell 
as it interacts with and dissipates into the surrounding material. Some cases are known in which 
the SNR has over-run a wind-swept shell, but the wind bubble invariably emerges irregular and 
fragmented in shape after its encounter with the supemova^^. The symmetric and well-defined 
appearance of our ring thus would imply that the supernova shock wave is still inside the ring. Our 
simulations show that this would require the supernova to be < 200 yr old, with bright radio and 
X-ray emission, contrary to observations. The lack of high-energy X-ray and non-thermal radio 
emission also excludes the possibility that the observed ring is produced by the emission of the 
blast wave from the progenitor's supernova explosion. 

The rejection of these various possibilities raises the issue of the energy source that is pow- 
ering the ring. We derived a limit for the required luminosity of the heating source by means of 
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radiation balance, using the measured temperature Tjr of the dust. If the luminosity is domi- 
nated by a source of temperature T*, then (see Supplementary Information for details) 

^ 2.7 X lO^^Tf^T-^R^ erg s"^ (1) 

where Tjr is measured in units of 100 K, is measured in units of 10^ K and R, the distance from 
the central heating source, is measured in units of parsecs. 

We note that the persistent X-ray luminosity of the SGR is only L2-iokeV ~ 2.0 — 3.5 x 
10^^ o?i5 erg <^ and, thus, cannot generate the observed infrared emission from the ring. 
The only viable heat source consistent with the observed properties of the ring appears to be ir- 
radiation from the stars in the nearby cluster. Figure la clearly shows that both the cluster and 
SGR 1900+14 are embedded in a diffuse, extended cloud of 24 /xm emission. Similar diffuse 
24 fjLm emission is also observed in association with the clusters containing SGR 1806—20 (ref.24) 
and Westerlund 1, implying that the extended emission is powered by the cluster stars. We note 
that, given the location and size of the 24 /xm nebula, the cluster might be more spatially extended 
than previously thought, with as-yet-undiscovered member stars. 

A possible mechanism for the creation of the ring is the formation of a dust-free cavity owing 
to the giant flare activity from SGR 1900+14, which would naturally explain why the ring is centred 
on the magnetar. The destruction of dust grains can occur either by sublimation due to the heating 
by ultraviolet^^ or X-ray^*' emission, or by grain charging due to the incident X-rays, which causes 
the dust grains to gradually shatter into smaller pieces until they are eventually destroy ed^^. The 
size and shape of the ring can constrain the energy output and degree of anisotropy of the flaring 
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event (see Supplementary Information for details). Using Equation (25) of ref.26 with the fiducial 
values and a = (where the flux at frequency v is proportional to i^~"), we find that the dust 
destruction radius matches well with the observed dimensions of the ring around SGR 1900+14 
for £■ > 6 X lO^^dfg erg. This estimate is close to the observed isotropic equivalent energy^^ in the 
initial spike of the 27 August 1998 flare; the latter, however, is only accurate for photon energies 
> 50 keV, whereas the estimated dust-destruction radius requires photons with energies of ~ 1 
keV. 

Alternatively, a previous, more energetic giant flare from SGR 1900+14, for example, similar 
to the 2004 December 27 giant flare from SGR 1806-20 (refs 27, 28), with E lO^^dj^ erg, could 
have carved the dust- free cavity during the magnetar's estimated spin-down age of ~ 1, 800 yr. For 
a rate of one flare per ~ 50 yr (ref.28), it is reasonable to expect at least one such event lO^^kyr yr 
ago, with tkyr ~ 1- Then, for the current location of the SGR still to be at the centre of the ring to a 
precision of 1", the proper motion of the magnetar (due to its birth kick velocity v±) and the proper 
motion of the dust-free cavity (with velocity Vc) should satisfy max(vc, v±) < 71di5t^^ km s~^. 

In light of the probable association between the SGR and the star cluster, we have re- 
examined the issue of the mismatch in extinction, which has been used as an argument against 
that connection^ ^ As our observations show, the cluster-SGR environment is characterized by a 
complex dust distribution on small spatial scales and dust destruction by the high-energy emis- 
sion of the SGR. Because the derived extinction values were based on measurements towards two 
distinct physical locations (the supergiants and the SGR) and are based on different tracers (gas 
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and dust) it is then not surprising that the resulting values differ. Similar differences have been 
observed for the Cas A SNR^^ and extragalactic gamma ray burst sources^^. In addition, we note 
that the visual absorption derived solely from the near-infrared photometry of the M supergiants, 
using both the published data^^ and new measurements extracted from the Two Micom All Sky 
Survey, is Ay = 10.4 — 14.4 mag, in comparison with Ay = 19.2 ± 1 mag, the value obtained 
when the optical I band measurement is included. Although this is puzzling and deserves further 
investigation, the various extinction measurements are clearly subject to large uncertainties and 
can potentially be reconciled with our conclusion that the SGR is a member of the star cluster. 
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Supplementary Material for "An infrared ring around the 

magnetar SGR 1900+14" 



Light echo scenario details 

If the SGR were surrounded by an extensive dust cloud, a substantial amount of the giant flare 
radiation will be converted into IR radiation. Given the extended nature of the dust cloud, light 
travel times across it must be considered. The two MIPS 24 //m images were obtained atti — 
7.10 yr and ta = 8.75 yr after the 1998 August 27 giant flare from SGR 1900+14. If the ring is 
from a light echo due to reprocessing of the giant flare photons then the time delay t is related to 
the location of the (assumed rapid) reprocessing by ct = r(l — cos 6'), where (r, 9, cj)) are polar 
coordinates with the SGR located at the origin and with the 2;-axis pointing to the observer. This 
locus of points from which reprocessed photons reach the observer simultaneously describes a 
paraboloidal surface centered on the location of the SGR and symmetric about our line of sight to 
the SGR, which grows linearly in time, in the sense that r/t is a function of the polar angle 9 alone. 
The lateral distance from the SGR to this surface (at 9 — 90°) is simply ct, which corresponds to 
cti «i 2.2 pc and ct2 ~ 2.7 pc for our two epochs of observation, where cti corresponds to the 
semi-major axis of the ring for d ~ 12.5 kpc. 

The main characteristics of such a light echo can be demonstrated by the simple and idealized 
case of a thin spherical shell of dust of radius R centered on the SGR. Its intersection with the 
paraboloidal surface described above is simply at cos^ = 1 — ct/R so that the apparent ring size 
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-^ring is 

R,,^,^Rsme = R^^^ (2- I) 0<|<2- (2) 

It initially grows super-luminally, its expansion slowing down until it reaches a maximal size of 
-Rring — Rati — R/c, and then begins to shrink increasingly fast until it disappears att — 2R/c. 
Its apparent velocity /Sapp in units of c is 



'7'/^ 0<|<2. (3) 

c dt let h- ^) ^ 



While /3app > 1 for much of the time, we also have P^pp <^ 1 at t a; R/c, i.e. very close to 
the time when the ring reaches its maximal size and begins to shrink. However, for two generic 
observation times one expects a fractional change of ~ (^2 — ^i)/^i in the ring size, which in our 
case corresponds to ~ 23%. Moreover, since the observed ring in our case is elongated with an 
axis ratio of ~ 2 : 1, the shell of reprocessing dust must clearly be non-spherical, and therefore the 
maximal size would also depend on the azimuthal angle (p, so that the apparent expansion velocity 
of the ring would never vanish everywhere at once, and one can expect fractional changes of the 
order of ~ (^2 — ^i)/^i ~ 23% in the ring shape or brightness between the two epochs, at least in 
some parts of the ring. Therefore, our strict limits on its exact location, and fractional chances in its 
size, shape, and brightness (of < 1 — 2%) are very hard to satisfy simultaneously in this scenario. 



Dust heating considerations 



The temperature of the dust can be used to determine the required luminosity of the heating source 
by means of radiation balance. If ex is the dust absorption efficiency, then its temperature is given 
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by a balance between the heat input due to the flux it absorbs and emits. If the flux it absorbs 
is dominated by a point source (which is also a useful approximation for a group of sources at 
distances larger than their separations) with luminosity L^, peaking at A*, at a distance R, then a 
(spherical, for simplicity) dust grain of radius a absorbs energy at a rate of 

-^absorbed = eA.TTa^- ^ , (4) 

and radiates away energy at a rate of 

-E'radiated = 47ra^eAij^aTj^ , (5) 
where Air ~ hc/kTi^. Equating the two rates gives: 

= 167ri?2 1^^^ aT/j, . (6) 

If the luminosity L^, is dominated by a source of temperature T^. (~ he/ k\^), and both e^. and cair 
are in the range where eA oc 1/A, then 

^A:, - -Qmi^M^y . (7) 



tx, T, ■ VlOOK/VlO^K 



and 



.40 / -'IR \ / \ ^ \ -1 



This basically sets a limit for the required heating source. Note that for a given source 
equation ^ implies that the dust temperature decreases rather slowly with the distance from the 
source, Tir oc R^"^/^ and could help explain why we do not detect a change along the ring in the 
apparent temperature of the ring emission. 
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Giant flare energetics and anisotropy 



The ~ 2 : 1 axis ratio of the ring suggests that if it was indeed formed by the initial spike of a 
giant flare, then the flare had an angular variation in its isotropic equivalent luminosity L or energy 
of a factor of > 4 (since the sublimation radius scales as L^/^ while the dust shattering radius 
scales as E^/"^, and the true axis ratio of the dust- free cavity could be somewhat larger than the 
observed axis ratio of the ring due to projection effects). Moreover, the relatively modest degree 
of anisotropy suggests that the true total energy output in the event that created the dust-free cavity 
was > 6 X lO^^d^g erg, which is comparable to the isotropic equivalent energy output for the 2004 
giant flare from SGR 1806-20. For the latter event it is hard to estimate the true energy output, 
since we have no observations that directly constrain the emission in directions other than our line 
of sight. 
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Figure 1 The infrared ring around SGR 1900+14. Spitzer Space Telescope imaging 
of SGR 1900+14 was acquired on 2005 September 20 using the InfraRed Array Camera 
(3.6-8.0 /xm), on 2005 Oct 1 1 witli tlie blue (16 ^im) peak-up imaging mode of the InfraRed 
Spectrograph (IRS), and on 2007 May 29 with the Multiband Imaging Photometer for 
Spitzer (MIPS; 24, TO/^m). Only the MIPS and IRS images are shown here. We assembled 
the basic calibrated data products into combined images using the Mosaicking and Point- 
source Extraction package. For the MIPS data, we first applied a differential flat field to 
correct for small (1-2%) instrumental artefacts. The final combined images (a-d) have 
effective exposure times of 31 s for IRS and 1 6 s for both MIPS channels. The white circle 
in the centre of each image marks the radio position of the SGR. For all images, north 
is up and east is to the left, a, MIPS 24 //m image scaled to show the diffuse emission 
enveloping both the SGR and the nearby cluster. The white circles labeled 'A' and 'B' 
indicate the positions of the red supergiants^^ located at the centre of the star cluster, 
b. Close-up of the MIPS 24 //m image shown in a, but scaled to highlight the ring-like 
structure of the extended emission, c, IRS 16/xm peak-up image covering the same field 
of view as b. d, MIPS 70 //m image showing emission associated with the enhanced 
areas along the minor axis of the ellipse also seen at 16 and 24/xm. The position of the 
ring is indicated by the black ellipsoid. The convolution of the 16 and 24 /xm images to the 
resolution of the 70 /im array indicates that the ring emission is too spatially confined and 
faint to be detectable at this wavelength. 
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